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ABSTRACT 

The shock-acceleration theory predicts a power-law energy spectrum in the test parti- 
cle approximation, and there are two ways to calculate a power-law index, Peacock's 
approximation and Vietri's formulation. In Peacock's approximation, it is assumed 
that particles cross a shock front many times and energy-gains for each step are fully 
uncorrelated. On the other hand, correlation of the distribution of an energy-gain 
factor for a particle is considered in Vietri's formulation. We examine how Peacock's 
approximation differs from Vietri's formulation. It is useful to know when we can 
use Peacock's approximation because Peacock's approximation is simple to derive the 
power-law index. In addition, we focus on how the variance of the energy-gain factor 
has an influence on the difference between Vietri's formulation and Peacock's ap- 
proximation. The effect of the variance has not been examined well until now. For 
demonstration, we consider two cases for the scattering in the upstream: the large- 
angle scattering (model A) and the regular deflection by large-scale magnetic helds 
(model B). Especially there is no correlation among the distribution of an energy-gain 
factor for every step in model A. In this model, we see the power-law index derived 
from Peacock's approximation differs from the one derived from Vietri's formulation 
when we consider the mildly-relativistic shock, and the variance of the energy-gain fac- 
tor affects this difference. We can use Peacock's approximation for a non-relativistic 
shock and a highly-relativistic shock because the effect of the variance is hidden. In 
model B, we see the difference of the power-law converging along the shock velocity. 

Key words: acceleration of particles - shock waves - methods:analytical - cosmic 
rays. 



1 INTRODUCTION 

Galactic cosmic rays with energies E ^ 10 15 eV, be- 
low the so called 'knee' in the cosmic ray spectrum, are 
thought to originate from shocks of supernova remnants. 
The shock-acceleration theory was proposed by various au- 
thors indfirjend*ntly_to_ex2l^m the o ri gin of Galactic c o smic 



triors independently to explain tnc origin ot Galactic cosmic 
ravs.lAxford, Leer fc Skadronl dl977f).lKrvmskvl (| l977l ).lBeiil 
|l97St ) and iBlandford fe Ostrikerl | l978l ') studied the shock- 



acceleration for non-relativistic shocks. Their theory relies 
on multiple scattering of charged particles which results in 
round-trips between the upstream and the downstream of a 
shock. In the test particle approximation, shock-acceleration 
theory predicts that a shock produces energetic particles 
whose energy spectrum obeys a power-law which explains 
the observation of the high energy cosmic rays well. A power- 
law is obtained by two ways. 1) Solving the convection- 
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diffusion transport equation with appropriate boundary con- 
ditions. 2) Considering an interaction of test particles with 
a shock front from the kinetic theory viewpoint. These in- 
dependent two theories predict the same result in spite of 
different approaches. 

Cosmic rays with energies E ~ 10 18 ' 5 — 10 20 ' 5 eV are 
generally believed to come from extragalactic origins, based 
on their harder spectrum, isotropic arrival directions on the 
sky and the fact that they are not confined by the Galactic 
magnetic fields. Some of authors c onsider cosmic r a y pro- 
duction for the relativistic shocks (jWaxmanl (| 19951 ) ; Ivietrl 
l|l995h ; lMirgrornl (Il995h). The relativistic shock-acceleration 
is first considered by fPeacockl (|l98lT ) considering relativis- 
tic kinematics based on th e Bell's theory taking ac co unt o f 
the beaming effect. Later, iKirk fc Schneider! i|l987T )l (Il988l ) 
solved the transport equation with appropriate collision op- 
erators. They focussed on the diffusion approximation. 

On the other hand, iKato fc Takaharal i|200lf) and lVietril 
formulated the shock-acceleration for the arbitrary 
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shock speeds without the diffusion approximation. Their 
theories inciude the non-diffusive effect which is important 
for reiativistic shocks and reproduce the same resuit in the 
non-relativistic limit. Both of them are constructed in two 
parts. First, they calculated P u (fiQ,^) (P ( j(/io,A t ) ) which 
denotes the conditional probability that a particle enter- 
ing into the upstream (downstream) along a direction /xo 
will leave toward the downstream (upstream) along a direc- 
tion jj,. They used some phenominological scattering models 
to calculate P u and Pd- Secon d, they calcula ted a power- 
law index using different ways. iPeacockl |l98ll ) calculated a 
power-law index assuming that particles cross a shock front 
many times and an energy-gain of each step is fully uncorre - 
lated (Peacock's approximation). iKato fc Takaharal d200lh 
used Peacock's approximation. On the other hand, IVietril 
solved the transport equation exactly and calculated 
a power- law index (Vietri's formulation). He considered cor- 
relation among energy-gains of a particle. He used his origi- 
nal formulation to calculate a power-law index. Vietri's for- 
mulation gives the different result from Peacock's approxi- 
mation in general but both appr oaches give the same result 
in the Newtonian limit (see also lBlasi fc Vietril ((2005)). 

To calculate a power-law index, we have to calculate 
P u an d Pd using some scattering models. IKato fc Takaharal 
l|200lf ) considered the large-angle scattering in both the up- 
stream and the downstream to calculate P u and Pd (Model 
A). The large-angle scattering model mimics the scattering 
in strongly turbulent fields. It is often supposed that the tur- 
bule nce can be st r ong in astrophysical shock environments 
(see lEllison et.al.1 fl990), and references therein). In this 
model, the initial information of the distribution function is 
lost when particles are scattered. So a n energy-gain factor of 
each s tep is uncorrelated. Moreover, iGallant fc Achterberd 
|l999I ) pointed out necessity to consider the regular deflec- 
tion as far as reiativistic shocks are concerned (Model B). It 
was understood that return of the particles to a shock sur- 
face from the upstream region can be warranted even in the 
absence of scattering, provided background magnetic fields 
are at an angle with the shock normal. This is due to the 
fact particles return to the shock surface from the upstream 
by the regular deflection before scattering occurs. 

In this paper, we examine how the variance of the 
energy-gain factor influences the difference between Vietri's 
formulation and Peacock's approximation. It is useful to 
know when we can use Peacock's approximation because 
Peacock's approximation is simple to derive a power-law in- 
dex. It is also important to understand Vietri's formula- 
tion and Peacock's approximation in detail because we can 
examine validity of the works in which Peacock's approxi- 
mation is used. First, we calculate P u and Pd in model A 
(where the large-angle scattering occurs at both sides of a 
shock) which is suitable for examining the effect of the vari- 
ance of the energy-gain factor because there is no effect of 
the correlation in this model. Then we calculate a power- 
law index using both Peacock's approximation and Vietri's 
formulation. Next, we consider model B where the large- 
angle scattering occurs in the downstream and a particle 
is deflected by large-scale magnetic fields in the upstream. 
We calculate probability functions and a power-law index 
using Peacock's approximation and Vietri's formulation in 
the same way as the case of model A. Finally, we examine 
how the power-law index derived from Peacock's approxima- 



tion differs from the one derived from Vietri's formulation 
to examine the effect of the variance. In model A, we see 
that the variance affects the difference between Peacock's 
approximation and Vietri's formulation, and explain what 
the effect of the variance means physically. To examine the 
above fact, we show the power-law index changing the shock 
velo city from a non-relativis t ic on e to a highly-relativistic 
one. iMorlino. Blasi fc Vietril (|2007l ) calculated a power-law 
index with 0.04 ^ T S [3 S ^ 10 considering the regular deflec- 
tion, where f3 s is the shock velocity and T s is the Lorentz 
factor of the shock velocity . We extend this calculation to 
the highly-relativistic range. 

The plan of this paper is as follows. In section 2 we 
briefly summarize the t heoretical framework introduced in 
IKato fc Takaharal l|200ll ). We also review Peacock's approx- 
imation and Vietri's formulation. In section 3, first, we con- 
sider the large angle scattering in both the upstream and the 
downstream. Next, we consider the regular deflection in the 
upstream and the large angle scattering in the downstream. 
We also discuss what makes Peacock's approximation inad- 
equate. Discussion and conclusion are presented in section 
4. 



2 METHOD OF CALCULATION 

In this paper we use the shock-acceleration formulation 
which is applicable to any value of the shock speed in the 
static situation and calculate a power-law index using Pea- 
cock'approximation and Vietri's formulation. We assume the 
test particle approximation and adopt the large-angle scat- 
tering. We also consider the regular deflection by large-scale 
magnetic fields, (we take the unit c = 1) 

2.1 Shock structure 

First, we have to determine a shock structure by solv- 
ing jump conditio ns and an equation of state, (e.g., 
iKirk fc Duffvl (|l999T )) Reiativistic jump conditions are writ- 
ten as 

TuPun u = FdPdTld, (1) 

r£/3„(e„ + p») = r 2 d /3 d (e d +p d ), (2) 

r„/3^(e u +Pu) +Pu = ^d0d{td +Pd) +Pd- (3) 

Number densities (n), pressure (p) and energy densities (e) 
are all measured in the comoving frame of the plasma we 
refer to, while the Lorentz factor (r) is measured in the 
shock frame. The indices 'u' and 'd' refer to the upstream 
and downstream plasmas respectively. 

For simplicity, we consider the case of a strong shock 
(i.e., upstream plasma is cold), so that p u = and e ~ n u m. 
We us e an equation of state which was introduced bv lSvngel 
(| 19571 ). 

The basic assumption is that the plasma consists of a 
single component with temperature T, 

e + p ^ G (5r)' < 4 > 

where G(x) — Ks(x)/K2(x) and K^,Kz are the modified 
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Bessel functions. If the downstream plasma can be well de- 
scribed as an ideal gas, Eq. (4) can be rewritten as 



, , nm 
e = pG\ 

P 



(5) 



Solving these equations, the solutions for (3d are obtained. 
(3 U and (3d are necessary for calculating P u and Pd- 



2.2 calculation of P u & P d 

In this subsection we briefly summarize the method of de- 
termining P u and P d in the case of both the large-angle 
scattering and the regular deflection. P u (Pd) is the condi- 
tional probability that a particle entering into the upstream 
(downstream) along a direction fio will leave toward the 
downstream (upstream) along a direction /j. 

First, we explain about th e case of the large-angl e 
scattering which was proposed bv lKato fc Takaharal (|200lT) . 
They assumed the energy of a particle measured in the fluid 
frame is conserved by scattering. The probability of displace- 
ment Az of the particle along shock normal is defined as 



p(Az,v, n)dAz — —e 
A 



~dAz. 



(6) 



Here, A denotes a mean free path of the particle measured 
in the shock rest frame, v is the speed of the particle and [i 
is the pitch angle cosine measured in the fluid frame. They 
use the scattering model that the pitch angle cosine of the 
particle [i is determined according to the probability den- 
sity function P a (fi,v). P u is independent of the initial pitch 
angle. This means the initial information is lost and there is 
no correlation between the initial distribution function and 
the last one. They consider the scattering is isotropic in the 
fluid frame: 



Pu = 



1 



Next, they defined the probability density function (p.d.f/ 
of the displacement for each step of the random walk as 

f{Az) = J P u p(Az,v,nW 

Jl u ^e-^d^(Az>0) 



/_7^e^d M (A*<0), 



(8) 



where u is the fluid velocity (e.g., Ud is the downstream fluid 
velocity in the shock rest frame). This p.d.f. expresses the 
probability of the particle which moves Az after the particle 
is scattered once. Then, we can calculate p.d.f. of scattering 
points at the m-th step, f m (Az). Finally, they introduced 
the density of scattering points summed over all steps that 
is written as 



n(Az) = J2 f m (Az 



(9) 



where n(Az) means that the particle is translated by Az 
before it crosses the shock front. By using the formulation 
mentioned above, we can estimate P u and P d after some 
detailed calculations. The results are: 



P u (/i , y) = Co(Ao) 



Pg\ 

A + A 



+ Ci(A ) 



P M A 
\ + L L 



(10) 



Pd(no,fJ.) = Cq(Ao) 



PaX 



+ Ci(A )P M A, (11) 



A + Ao 

where Co, Ci, Co, C{ are the functions which depend on the 
shock speed and the mean free path. Lp is t he dif fusion 
length which is defined by IKato fc Takaharal i|200ll ). The 
total return probability in the downstream, Pr, is not unity, 
and given as 



Pr(vo) = 



Pd(/J.o, n)dfi. 



(12) 



They also calculate the flux across the shock front. The re- 
sults are 



</)du(n) = 



<t>ud{n) 



Pu\ 



(13) 



(14) 



g+(L D ) A + Ld ' 

as fo r the estimation of h- and g+ , see IKato fc Takaharal 
2001)). 4>d u {(pud) is the flux entering to the upstream 
(downstream). 

Next, we explain about the case of the regular deflec- 
tion in the upstream. The regular deflection is important 
for a relativistic shock in the upstream. In the case of the 
regular deflection, we have to solve the equation of motion 
in la rge-scale magneti c fields l|Gallant fc Achterberd (Il999l ) 
and lAchterberg et.al.1 (|200lh 'l . Here, we choose the shock 
normal in the z-direction and assume the magnetic fields 
is parallel to the x-direction. We calculate the pitch-angles 
jtii when particles cross the shock front from the upstream. 
First, we decide the initial azimuthal/pitch angles fj,o and 
0o- The particle's velocity is decided as 

(3 X = cos 4> sin 0, fjy — sin sin 6, /3 Z — cos 0. (15) 

Next, we solve the equation of motion which is written as 



d/3 
~dt 



O s (/3 x b) 



(16) 



where q — Ze is charge of the particle, E is energy and Q g — 
ZeBc/E is a gyration frequency in magnetic fields B—B b. 
We solved this equation numerically using the above velocity 
as the initial condition. The upstream residence time t u is 
obtained by considering z(t u ) = z 3 — [3 s t u . Here, z(t u ) is the 
position of the particle and z s is the position of the shock 
front. Finally, we calculate the pitch-angle /ii using Eq. (15) 
again, /ii is determined uniquely as a function of fio- We can 
express P u using four angles as 



P u {po,y) = (2tt) 1 5(p - Hi((j,o, </>o)) 
x5(cj> - (j>i(p,o, <t>o))- 



(17) 



2.3 Peacock's approximation 

iPeacock] |l98li) developed the calculating method to ob- 
tain a power-law index which was originally used by iBelll 
(1978). Peacock's approximation can be applied to a rela- 
tivistic shock. He considered iVo particles crossing from the 
upstream to the downstream with an initial energy Eq. We 
use a return probability Pr((j,o) that a particle will even- 
tually return to the upstream. The particle's energy is in- 
creased by a energy-gain factor G(/j,o,(J,), where /io and p 
are the initial pitch angle and the last pitch angle for one 
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step. After k cycles, number of particles that remains in the 
upstream (N) is expressed as 



where 



(Ph) 



PjlOio)0«d(jUo)<fyio. 



(18) 



(19) 



Here, (Pr) is the return probability averaged over fio with 
the weight of the flux crossing the shock with various values 
of fio. As particles cross and re-cross the shock, the distri- 
bution of energies is broadened. Peacock's approximation 
assumes the numbers of cycles are large enough to use the 
central limit theorem. It is also necessary that the energy- 
gain of each step is fully uncorrelated. If particles cross the 
shock front many times, the distribution of energies can be 
expressed by a Gaussian by a central limit theorem. In such 
a situation, we can approximate particle's energy is ampli- 
fied at the same rate per one cycle and the effect of the 
variance is neglected. We can calculate the power-law index 
only using the averaged energy-gain factor. After k cycles, 
the particle energy is given by 



In 



!■)=/.. in,; 



(20) 



where G is the energy-gain factor and expressed as 
G((i ou t,n) = ( 1 ~ s ^y°^ t ^ ■ The averaged InG is defined as: 



(In G) 



1 r-u.d 

dfi / dfi ou t 

■u d J -1 



•' '« ! /»..., ' - /' ) In ( ^ i Tyjf ) ' ' /! '"" ; 1 



dfl ut<\>du {[lout) 



(21) 



Here, V r 



indicates the relative velocity of the 



upstream fluid with respect to the downstream fluid. From 
Eq. (18) and (20), the integrated energy spectrum is written 
as 



1 (*L 

n \N J (InG) ~\E 0j 
Thus the differential energy spectrum is obtained as 



where 



f(E)d i p oc E~ s dfp, 

MP*) 



s = 3 



(In G) 



(22) 



(23) 



(24) 



2.4 Vietri's formulation 



Ne xt, we e xplain derivation of a power-law spectrum given 
bv IVietril J2003). He derived the relativistically covariant 
equation for the distribution function of particles accel- 
erated at a shock, which is applicable to a relativistic 
shock. His formulation gives the exact power-law index and 
does not assume unc orrelation among various energy-gains. 
iBlasi fc Vietrll (|2005h solved the transport equation provid- 
ing the boundary condition for the flux which crosses the 
shock front to the downstream (4> u d) (upstream (cfidu)) by 




1000 



Figure 1. Return probability vs the shock speed for model A 
(solid line) and model B (dotted line). We assume large-scale 
magnetic fields are present in the upstream, with a direction per- 
pendicular to the shock normal when we consider the regular 
deflection (model B). 



using particles which crosses the shock front to the upstream 
{<j>du) (downstream (4>ud)), and they derived P u and Pd- In 
this paper, we use P u and Pd which are calculated in section 
2.2. The power-law index is given by this boundary condi- 
tion: 



4>ud{y) = 



d/loutPuifJ'out, fl) 



1 - Vrfl 
1 — V r flout 



4>du{flout 



(25) 



This equation can be integrated over the whole range of 
/j, and divided by the whole flux entering into the upstream, 
which gives 



J_ 1 d dfl ou t4>du{Hout) 



dfi 



dfiout 



Pu(flout, fl) 



1 - Vrfl 
1 — Vrflout 



4>du(flout) 



dfl out 4>du (flo u t 



(26) 



The left term is the inverse of the averaged return probabil- 
ity from the downstream. The right term is the average of 
the (s — 3) power of the energy-gain factor G. Eq. (26) can 
be rewritten as 



(P fl )(G s " 3 ) = l. 



(27) 



In Newtonian limit, we can approximate (G s 3 ) ~ 
{G) s - 3 because G- 1 « 1 for a non-relativistic shock. Thus 
the energy spectral index is given by the similar form as 
Eq. (24). Vietri insisted that (G s_3 ) can not be approxi- 
mated by (G) s_3 in general for a relativistic shock, contrary 
to the Peacock's argument. 
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Figure 2. The averaged energy-gain factor for model A (top 
panel) and model B (bottom panel). We can see the averaged 
energy-gain factor of model A is much larger than that of model 
B. 



3 RESULTS 

In this section we calculate a power-law index using Pea- 
cock's approximation and Vietri's formulation. Then we ex- 
amine the characteristic of the power-law index and how 
the power-law index derived from Peacock's approximation 
is different from the one derived from Vietri's method to 
examine the effect of the variance. We consider model A 
(where the large-angle scattering occurs in both the up- 
stream and the downstream) , and model B (where the large- 
angle scattering occurs in the downstream and particles are 
deflected by large-scale magnetic fields in the upstream). 
We need the averaged chance probability for a particle to 
return from the downstream to the upstream per a crossing 
cycle (Pr) (Eq. (19)) and the averaged energy-gain factor 
(InG) (Eq. (21)) to calculate the power-law index by Pea- 
cock's approximation. In Fig. 1 we show the averaged return 
probability. In the both cases of model A and model B, the 
return probability converges along the shock velocity. This 
is because the shock velocity converges to one third of the 
speed of light in the down stream when we use the Synge 
equation as the equation of the state. So the probability 
for particles to catch up with the shock front converges. In 
Fig. 2 we show the averaged energy gain factor. The value 
increases monotonically in model A (Fig. 2 top). On the 
other hand, the value converges in model B (Fig. 2 bot- 
tom). Particles get more energy by the shock-acceleration 
when the shock velocity becomes faster and the scattering 
angle becomes larger. In model A, particles are scattered at 
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Figure 3. The power-law index vs. the shock speed. The power- 
law indices in model A and model B are calculated by Peacock's 
approximation and Vietri's formulation. The above two lines cor- 
respond to model B, while bottom two lines correspond to model 
A. 



a large angle and get large energy as the shock moves with 
the large velocity. In model B, the deflection angle becomes 
smaller as the shock velocity becomes large. Both of the two 
effects (a large velocity and a small scattering angle) cancel 
out each other and the energy amplification converges. In 
Fig. 3, we show the power-law index. The power-law index 
becomes harder as a shock moves fast in model A. This is 
because the averaged energy-gain factor increases monoton- 
ically. On the other hand, the power-law index converges in 
model B. This is due to the convergence of the averaged re- 
turn probability and the averaged en ergy-gain factor. This 
conve rgence is the same behavior as iBednarz fc Ostrowskil 
(1998). They studied about the case of the small-angle scat- 
tering and the power-law index converges as the shock ve- 
locity becomes large. The power-law index converges to 4 
for a non-relativistic case which is consistent with the result 
of the diffusive shock-acceleration. 

Next, we consider the difference between the power- 
law index derived from Peacock's approximation and Vi- 
etri's formulation. The power-law index converges in model 
B, but the converged values are different between Peacock's 
approximation and Vietri's formulation. So Peacock's ap- 
proximation is not suitable for relativistic shocks. On the 
other hand, in the case of model A, the difference of the 
power-law index is the largest when T s (3 3 is about 3, and the 
index becomes nearly equal when T s f3 s is greater than 10. 
We can use Peacock's approximation for highly-relativistic 
shocks. We also show that the power-law index calculated 
by Vietri's formulation is harder than the one calculated 
by Peacock's approximation in both cases of model A and 
model B. We explain why the power-law spectrum becomes 
hard in the discussion. 

Next we see what makes the difference between Pea- 
cock's approximation and Vietri's formulation. First, trans- 
form Eq. (24) and Eq. (27) as 



ln(P r (/3 sl ))+ln<G sl - J (/? sl ))=0, 
HPrO&J) + («a - S)(\nG(p si )) = 0, 



(28) 
(29) 



where Si(s2) is the power-law index derived from Vietri's 
formulation (Peacock's approximation) and /3 S1 is the shock 
velocity which satisfies Eq. (27). Peacock's approximation 
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Figure 4. The encrgy-gain factor's distribution function P(G) 
in model A (top panel and middle panel) and model B (bottom 
panel). The results are plotted in the cases of several values of 
the shock velocity. In model A, the plotted distribution function 
is rescaled (multiplied by 10) both for (3 S = 0.3 and f3 s = 0.5. At 
the low velocity (top panel), the distribution function does not 
looks like a rectangle. At the high velocity (middle panel), the 
distribution function looks like a rectangle and we use the rect- 
angular distribution function. In model B, we plotted the distri- 
bution function when <j> is ^ir. The plotted distribution is rescaled 
(multiplied by 10) both for j3 s = 0.5 and j3 s = 0.99. In this case, 
the distribution does not look like a rectangle and we can not use 
approximation. 



(Eq. (24)) and Vietri's formulation (Eq. (27)) coincide with 
each other as long as the following equation is satisfied: 

ln<G sl - 3 (/3 sl )> = ( Sl -3)(lnG(/3 sl )>. (30) 

Now we introduce the value defined as 

D(J3 tl ) = ln(G sl - 3 (/3 sl )> - (si - 3)(lnG(/3.i)>. (31) 

We use D((3 S1 ) as the function which indicates the difference 
between Peacock's approximation and Vietri's formulation. 
The averaged value is calculated using the energy-gain factor 
distribution as follows, 



(G) = J GP(G)dG, 



(32) 
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Figure 5. D and D app vs. the standard deviation. The approxi- 
mation gives the contour with the power-law index fixed. It could 
be understood that the effect of the variance becomes weak as 
the power-law index approaches 3. 



where 



P(G) 



r u 1 - 

J -UA \ 



Vrfl 



(33) 



D becomes if Peacock's approximation and Vietri's for- 
mulation give the same result. Particularly, D approaches 
when s approaches 3. D becomes large if the difference 
between Peacock's approximation and Vietri's formulation 
becomes conspicuous. 

Now let us see the effect of the variance in D. How- 
ever, it is difficult to know the relation between D and the 
variance. Here, we use the simple model of the distribution 
function which shows the distribution function (Eq. (33)) 
approximately. In Fig. 4 the energy-gain factor's distribu- 
tion function P(G) is plotted for the various shock velocity. 
In the case of model A, the distribution becomes wider and 
rectangular as the shock moves fast. Then we use the rect- 
angular distribution function for highly-relativistic shocks 
in model A. The rectangular distribution function is written 
as, 



Papp{G^j — 



1 

L-l 

other 



1 «: G < L 



(34) 



On the other hand, the distribution does not become rectan- 
gular and seems unsuitable to apply this simple distribution 
function in model B. L>, which is calculated using the simple 
distribution model, is written as, 



D a pp(L,si) = In 



1 



si - 2 
si — 3 



L 



si —2 



L-l L-l 
(LlnL — L + 1) . 



(35) 



We expect that the variance which is calculated by the 
true distribution function (Eq. (33)) corresponds to the one 
which is calculated by the rectangular distribution. L has a 
relation to the variance a 2 as a = (L- l)/2\/3. Using this 
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Figure 6. The variance of the energy-gain factor's distribution 
function vs. the shock speed for model A (top panel) and model 
B (bottom panel). 



relation, D app {L,S\) is rewritten as, 



In 



si - 2 
si — 3 



(2Vga + 1) S1 ~ 2 _ 1 
(2V3a + 1) \n(2V3a + 1) - 2VSa 



(36) 



D ap p has the convergence value when a becomes infinity. 
The convergence value is written as, 

1 



D app (oo, si) = In • 



si 



si - 3. 



(37) 



We examine how D depends on the variance in model 
A. In Fig. 5, we can plot D and D app as the functions of the 
standard deviation (i.e., the variance) instead of (3 3 i because 
the variable increases monotonically as /3 s i increases. We 
show the relation between the shock velocity and the vari- 
ance in Fig. 6. There is a one-to-one correspondence between 
the variance and the shock speed, a and s is written as the 
functions of one variable as, a = <r(r/3), s = s(T(3) — s(a). 
We calculate D using si which is calculated by Eq. (27). We 
plot the contour of D app fixing si mathematically in order 
to examine the effect of the variance. D app decreases as si 
approaches 3 for the fixed variance. This means the effect of 
the variance becomes weak as s\ approaches 3. D apP gives 
inadequate approximation when the shock speed is small. In 
fact, D app caluculated by s\ — 3.53 gives the same value as 
D although D app calculated by si —4 should give the same 
value as D. For highly-ralativistic, shocks, the approxima- 
tion becomes good. D changes along the contour which is 
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Figure 7. D is plotted both in model A and model B. We can 
interpret that D is derived from the effect of the variance in model 
B. In model A, D has maximum value . In model B, D converges 
as the shock speed becomes large. 



derived from the rectangular approximation at first because 
the power-law index does not change when the shock speed 
is low. D crosses the contour as the shock moves fast because 
the power-law index approaches 3. This behavior shows that 
the effect of the variance becomes weak by the approach of 
s to 3. 

In Fig. 7, D is plotted both in model A and model B 
as the functions of F S /3 S . In model A, D becomes the largest 
when F s /3 a is about 3, and the variance is effective. Although 
the variance increases monotonically, D decreases for highly- 
relativistic shocks. This is because the effect of the variance 
becomes weak as the power-law index approaches 3. D be- 
comes for non-relativistic shocks because the variance be- 
comes small as we can see from Fig. 6. In model B, we can 
interpret behavior of D as explained before with respect to 
the variance and the power-law index. D increases as the 
variance increases at first because the variance increases as 
the shock velocity becomes large (see Fig. 6) . After the vari- 
ance converges, D continues to increase until F S /3 S becomes 
as high as 30. This is because the power-law index increases 
and the effect of the variance becomes strong until F a f3 s 
becomes as high as 30. For non-relativistic shocks, D also 
bec omes in model B . 

iKato fc Takaharal i|2005t) considered the effect of the 
variance in their shock-acceleration formulation. The effect 
of the variance also becomes negligible for a non-relativistic 
shock in the case of their formulation. 



4 DISCUSSION & CONCLUSION 

In this paper, we examine how the power-law index derived 
from Peacock's approximation is different from t he one de- 
rived from Vietri's formulation. iPeacockl (Il981[) gives the 
approximate power-la w index and IVietril (|2003l ) gives the 
exact power-law index. IVietril (|2003l ) claimed that Peacock's 
approximation and Vietri's formulation give the different re- 
sult in ge neral, but they did no t explain the reason in detail 
(see also iBlasi fc Vietril l|2005l )). The effect of the variance 
has not been studied well until now and we have examined 
the effect of the variance. 

First, we conclude the difference of the power-law in- 
dex derived from Peacock's approximation and Vietri's for- 



8 Junichi Aoi, Kohta Murase and Shigehiro Nagataki 



mulation. Wc considered two cases for the scattering in the 
upstream to model multiple shock crossing: the large-angle 
scattering (model A) and the regular deflection by large- 
scale magnetic fields (model B). We consider only the large- 
angle scattering in the downstream. In model A, the power- 
law index derived from Peacock's approximation is different 
from the one derived from Vietri's formulation for a mildly- 
relativistic shock. The difference decreases as the shock ve- 
locity approaches the highly-relativistic velocity. In model 
B, the difference becomes larger along the shock velocity at 
first. Finally, the difference converges when the shock ve- 
locity is mildly-relativistic. From the above result, we can 
use Peacock's approximation for a highly-relativistic shock 
when we consider model A and the difference between Pea- 
cock's approximation and Vietri's formulation is important 
when we consider model B for a relativistic shock. More- 
over, we see the reason why Peacock's approximation and 
Vietri's formulation give the different results. We conclude 
the variance in the distribution of the energy gain factor 
affects the difference between Peacock's approximation and 
Vietri's formulation in model A. 

In Fig. 3, we show that the power-law index calculated 
by Vietri's formulation is harder than the one calculated 
by Peacock's approximation in both cases of model A and 
model B. This tendency can be understood as follows: Let us 
consider the number of particles with Eh, which have expe- 
rienced k cycles in an average (see Eq. (20)). If there is the 
variance of energy-gain factor's distribution, there should 
be contribution to the number from those that have expe- 
rienced less than fc-cycles and more than fc-cycles. It is ap- 
parent that the contribution from the former is larger than 
the latter, because the energy spectrum obeys a power-law 
with index larger than 3. As a result, the spectrum becomes 
hard when the effect of the distribution width is taken into 
account. As above, we can understand why the power-law 
index becomes hard as the variance becomes large. However, 
this interpretation can not be suitable when the variance is 
larger than some value. The effect of the variance converges 
as we show in Fig. 5. What's more, we can interpret that 
the effect of the distribution's variance becomes weak as the 
power-law index approaches 3 because the fraction of par- 
ticles that are accelerated by the mean energy-gain factor 
increases. 

Next, let us explain the general characteristics of the 
power-law index when we change the scattering model. We 
found the power-law index converges as the shock veloc- 
ity increases in model B. Th is is the same characteristic as 
iBednarz fc Ostrowskil (I1998T ) who studied about the case of 
the small-angle scattering, although the spectrum derived 
from the large-angle scattering is steeper than that of the 
small-angle scattering. On the other hand, the index de- 
creases monotonically in model A. For the limit of non- 
relativistic shocks, the result is consistent with the diffusive 
shock-acceleration theory. 

The difference of the power-law index between model 
A and model B is large for a relativistic shock. The differ- 
ence is small if one uses the small-angle scattering model 
instea d of the large-angle scattering (e.g.. [Achterberg et.al.1 
This is because the averaged energy-gain factor is 
very large in model A. The energy-gain factor can be writ- 
ten as G(^ ut,fJ.) = ( ^-v^"' ) • Thus, the energy-gain fac- 



tor can be large for a highly-relativistic shock in the case of 
the large-angle scattering since the difference between fi out 
and is large in the case of the large-angle scattering. On 
the other hand, the energy-gain factor remains the order of 
unity even for a highly-relativistic shock both in the case of 
the small-angle scattering and the regular deflection. 

Peacock's approximation and Vietri's formulation give 
the different results in general. To see how Peacock's ap- 
proximation is good, of course it is also important to inves- 
tigate how many times scatterings should occur enough to 
use the central limit theorem and whether there is correla- 
tion among energy-gains for every step. Such consideration 
about these effects should be done not only for relativistic 
shocks but also for non-relativistic shocks. For example, we 
probably can not use a central limit theorem even for non- 
relativistic shocks in low energy range. This consideration 
might be useful to discuss the boundary between the supra- 
thermal particle and the non-thermal particle. The variance 
might also be considered to calculate the power-law index 
as we show in this paper. 

We have considered the case of the large-angle scatter- 
ing in this paper. We will consider the case of the small-angle 
scattering in the future work. 
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